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Abstract. Following the works by Wiegmann-Zabrodin, Elbau-Felder, Hedenmalm-Maka- 
rov, and others, we consider the normal matrix model with an arbitrary potential function, 

■ and explain how the problem of finding the support domain for the asymptotic eigenvalue 
I density of such matrices (when the size of the matrices goes to infinity) is related to the 

problem of Hcle-Shaw flows on curved surfaces, considered by Entov and the first author in 
1990-s. In the case when the potential function is the sum of a rotationally invariant function 

■ and the real part of a polynomial of the complex coordinate, we use this relation and the 
conformal mapping method developed by Entov and the first author to find the shape of 
the support domain explicitly (up to finitely many undetermined parameters, which are to 
be found from a finite system of equations). In the case when the rotationally invariant 

■ function is (3\z\'^, this is done by Wiegmann-Zabrodin and Elbau-Felder. We apply our 
I results to the generalized normal matrix model, which deals with random block matrices 

■ that give rise to *-representations of the deformed preprojective algebra of the affine quiver 
of type Am-i- We show that this model is equivalent to the usual normal matrix model in 
the large N limit. Thus the conformal mapping method can be applied to find explicitly 

[ the support domain for the generalized normal matrix model. 

, Key words: Hele-Shaw flow; equilibrium measure; random normal matrices 
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. The normal matrix model became a focus of attention for many mathematical physicists after 

the recent discovery (see e.g. [HI EJ [TJ [8]) of its unexpected connections to the 2-dimensional 
dispersionless Toda hierarchy and the Laplacian growth model (which is an exactly solvable 
model describing free boundary fluid flows in a Hele-Shaw cell or porous medium). The original 
normal matrix model contained a potential function whose Laplacian is a positive constant, but 
later in [12], Wiegmann and Zabrodin considered a more general model, where the potential 
function was arbitrary. This is the model we will consider in this paper. 

In the normal matrix model with an arbitrary potential function, one considers the ran- 
dom normal matrices of some size with spectrum restricted to a compact domain and 
probability measure 

PN{M)dM = Z^^exp(-iVtrI^(M))dM, 



1 Introduction 
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where dM is the measure on the space of normal matrices induced by the Euchdean metric 
on all complex matrices, Vl^ is a potential function (a real function on C with some regularity 
properties, e.g. continuous), and Zjy is a normalizing factor. 

In the original works on the normal matrix model, the potential was 



where P is a complex polynomial of some degree d, and f3 a positive real number. For this type 
of potential, it was shown in the works [HI [H [71 [8] (and then proved rigorously in [3]) that under 
some conditions on the potential, the asymptotic density of eigenvalues is uniform with support 
in the interior domain of a closed smooth curve. This curve is a solution of an inverse moment 
problem, appearing in the theory of Hele-Shaw flows with a free boundary. Thus, applying the 
conformal mapping method (see [10] and references therein), one discovers that the conformal 
map of the unit disk onto the outside of this curve which maps to c« is a Laurent polynomial 
of degree d. This allows one to find the curve explicitly up to finitely many parameters, which 
can be found from a finite system of algebraic equations. 

In [12], Wiegmann and Zabrodin generalized this analysis to an arbitrary potential function. 
They showed that the density of eigenvalues is the Laplacian of the potential function, and the 
eigenvalues are concentrated in the domain which can be determined from an appropriate inverse 
moment problem. This was proved rigorously in the paper [5], which extends the Elbau-Felder 
work to the case of an arbitrary potential. 

One of the goals of the present paper is to use the generalized conformal mapping method, 
developed in [4] by Entov and the first author for studying Hele-Shaw flows with moving bound- 
ary for curved surfaces, to calculate the boundary of the region of eigenvalues explicitly in the 
case when 



where $ is a function of one variable. In this case, the conformal map of the disk onto the 
outside of the curve is no longer algebraic, but one can still give an explicit answer in terms of 
a contour integral. 

Another goal is to extend the above results to the case of generalized normal matrix model. In 
this model, we consider block complex matrices of a certain kind with commutatation relations 
similar to the definition of a normal matrix; they give rise to *-representations of the deformed 
preprojective algebra of the affine quiver of type Am-i- We prove that the problem of computing 
the asymptotic eigenvalue distribution for this model, as the size of the matrices goes to infinity, 
is equivalent to the same problem for the usual normal matrix model. This allows one to find 
the boundary of the eigenvalue region explicitly if the potential is given by ([T]). 

The structure of this paper is as follows. In Section [2l we state some basic facts about 
the normal matrix model. In Section [3l we define the generalized normal matrix model, and 
write down the probability measure in this model. In Section [U we recall some facts about 
the equilibrium measure and explain that the asymptotic eigenvalue distribution tends to the 
equilibrium measure in the normal matrix model and the generalized normal matrix model. In 
Section [5l we use the singular point method from [U to reconstruct the boundary of the 
support domain of the equilibrium measure. 

2 The normal matrix model with an arbitrary potential function 

Let D be a domain in the complex plane C. We consider the set 



W{z) =/3|z|2-P(z)-P(z), 




(1) 



Af{D) = {M G Mativ(C) I [M,Mt] = 0, spectrum(M) C D} 
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of normal matrices with spectrum in the domain D. Let dM be the measure on J\f{D) induced 
by the Euclidean metric on Mat7v(C). It is well known (see e.g. OH]), that in terms of the 
eigenvalues this measure on AA(C) is given by the formula 

N 

dM = JJ - I ^ n d^-^idC/, 

l<i<j<Ar i=l 

where M = [/diag(2;i, . . . ,Z]\])U\ U G U{N), and dU denotes the normalized [/(A^)-invariant 
measure on the flag manifold U{N)/U{1)^ . 

Now let : C ^ M be a continuous function. If M is a normal matrix, then we can define 
W{M) to be diag(l^(2;i), . . . , W{zi\f)) in an orthonormal basis in which M = diag(zi, . . . , zn). 
It follows from the above that the probability measure on M{D) with potential function W is 
given by 

N 

Pjv(M)dM = Z-ie-^E>^(^') n \zi-z,?\{d^ZidU, (2) 

l<i<j<N i=l 
N 

where Zn = J^jv e^^^''^^^'^ W \zi — Zj\'^ W d^Zi. Here we assume that the integral is 

l<i<j<N 1=1 

convergent (this is the case, for instance, if D is compact). 



3 The generalized normal matrix model 

3.1 Generalized normal matrices 

Let us consider the following generalization of normal matrices. Let m > 1 be an integer. For 
a fixed collection A = (Ai, . . . , Am) of real numbers such that Yli Aj = 0, and a domain D, we 
define MmiX, D) to be the subset of A G MatmAr(C) satisfying the following conditions: 
for any A G Nm{\, D), 

• If Aij, i, j = 1, . . . ,171 are N x N blocks of A, then Aij = unless j — i = 1 mod m; 

• The spectrum of A12A23 • ■ ■ Ami is contained in D; 
XiIn \ 

, where In is the identity matrix of size A^. 



[A, At] 



Am-^TV / 



Note that J\fi{0,D) = J\f{D), thus elements of A/'m(A,-D) are a generalization of normal 
matrices. We will thus call them generalized normal matrices. 

Remark 1. Generalized normal matrices are related in the following way to quiver representa- 
tions. Let Q be the cyclic quiver of type Am-i, and Q its double. Let nQ(A) be the deformed 
preprojective algebra of Q with parameters A (see [2]). By definition, this algebra is the quotient 
of the path algebra of Q by the relation ^ [a, a*] = ^ AjCj, where are the vertex idempotents. 

aSQ 

The algebra IIq has a *-structure, preserving and sending a to a* and a* to a. It is easy 
to see that Mm{X,D) is the set of all matrix *-representations of IIq of dimension N5 (where 
(5 = (1,1,...,!) is the basic imaginary root) such that the spectrum of the monodromy operator 
ai • • • am is in D. 
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Denote Au+i by Ai. The group [/(A^)™ = U{N) x U{N) x • • • x U{N) acts naturally on 



Mm{X, C) by the formula 



('S'l ) • • • 1 Sm)A 



SiAiS. 



S2A2SI 



\ SmAmS\ 



Sm—l-^m—lS; 



J 



We have the following lemma, which is a generalization of the fact that a normal matrix 
diagonalizes in an orthonormal basis: 

Lemma 1. For any element A G Mmi^^D), we can find an element {Si, . . . , Sm) £ U{N)'^ 
such that 



I 



(S'l , . . . , Sm)A 



\ 



D2 



n 



m—l 



J 



where Di are diagonal matrices. 

Proof. From the definition, we have — = XiIn, where the index is considered 

m 

modulo m, and ^ Aj = 0. 

i=l 

Now consider a collection of A^-dimensional unitary spaces {Fj}"!]^, and let us regard Ai as 
a linear map Ai : T^+i ^ V^. So AiAj is a hermitian endomorphism of Vi. 

Now suppose that v is an eigenvector of AiAl with eigenvalue u. We claim that then Ai-iv 
(if it is nonzero) is an eigenvector of Ai-iAl_-^ with eigenvalue u — Aj. Indeed, 

Ai^iAl_-^Ai_iv = Ai_i{AiAl - X^)v = {v - A^Ai-iw. 

Thus, denoting by Vi{v) the eigenspace of AiA\ with eigenvalue we find that Ai-i : Vi{u) — > 
Vi-i{i' — Xi). Since Vi = (BueRViii^) (as AiAj is hermitian), it suffices to prove the lemma in the 
case when AiA\ is a scalar in V^, in which case the statement is easy. ■ 

3.2 The Euclidean measure on generalized normal matrices 

First, let us consider the = 1 case. Pick real numbers ai such that Xi = ai — ai-i, and let 

m 

Q{x) = n + ^i)- ^ matrix 

i=l 



( 



A 



r2e 



1^2 



iQ-n 



\ 



(where rj > 0, 9j G [0, 27r)) is in Afm{X, C) if and only if 
'^i - ^m. = Ai = ai - am, 
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r| - ri = A2 = 02 - ai, 



A A \ 



Thus to each A G Mm{\, C) we can attach a real number x 



ai, which is independent 



of i, and a complex number = H fj^^^- It is easy to see that the point {z, x) belongs to the 
surface 



S = |(z,x) G C X ]R|x + ai > Vi, zz = JJ(x + ai)}. 



1=1 



Moreover, it is clear that any point of S corresponds to some A, and two matrices A, A' 
giving rise to the same point (z, x) are conjugate. This implies that we have a bijection between 
the equivalence classes in J\fm{X, C) under the action of U (I)'" and points of S. Writing z = re*^, 
we see that x, 6 are coordinates on S, so we may write the Euclidean measure on Mmi^i'C) 
using the coordinates x, 9. 

Theorem 1. The Euclidean measure on Mm{\,C) for N = \ is: 



l_t = -^Q'{x)dxdedU, 
where dU is the Haar measure on 



Proof. We have = i/x + Oj. Thus the parametrized real curve {(ri(x),... ,rm,(x))|x G 
has length element 



ds = rS^ r'Ax^dx = - , dx. 



This implies that the Euclidean measure on J\fm.{^, C) is 



2 \l ^-^ X + Oi 
V * 

^y/Q'{x)dxd9i ■ ■ ■ dOrn = ^VQ'{x)dxd9dU, 



as desired. 



■d9r^ 



Let us now consider the case of general N. From Lemma [H we know that under the action 
of U{N)"^, the equivalence class of A G A/'m(A, C) can be represented by m diagonal matrices 



M 



where Di = diag(2;j^, . . . , zj^). 
From the definition 



DiDI - DlDrn = All 



N, 



\ 



m—l 



(3) 
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D2DI - DlDi = X2IN, 
So we have ZjZj — Zj'^^Zj'^^ = cti — di+i- Let xj = ZjZj — ai and Zj = Hi then we have 

i 

Thus ((zi, xi), . . . , {z]\f, X]\f)) is a point on T,^ /Sn- 

Similarly to = 1 case, it is easy to show that this gives rise to a bijection between conjugacy 
classes of elements of A/'m(A, C) and points of S^/S'tv. Using this fact and combining the method 
of computation for usual normal matrices with the = 1 case, one gets the following result. 

Theorem 2. The Euclidean measure on A/'rrt(A, C) has the form 
dM = ^ JJ VQ'(xi) JJ \zi - Zjf dxi • • • dxivd^i • • • dONdU, 

i i<j 

where dU is the normalized invariant measure on U{N)"'' /U{\)^ . 

Proof. At first, consider the subset A/'m'^^(A,C) of A/'m(A,C) consisting of the elements M of 
the form ([3]). Then by Theorem[Tl the measure on Afm^^i^, C) induced by the Euclidean metric 
is the product measure: 

/^diag = ^ H VQ'i^i)'^^! ■ ■ ■ dx^dOi ■ ■ ■ d^AfdC/diag, (4) 

i 

where df/diag is the Haar measure on U{1)^"^ /U{1)^ . 

Now consider the contribution of the off-diagonal part. Consider the elements 

{vij = Eij - Ej^i, Wij = \^{Eij + Ej^i) \ ^i < j ^ N} 

of the Lie algebra of U{N). 

Let Vi,j,k, ^i,j,k be the derivatives of {exp{tVij))kM and {exp{tWij))kM at t = 0, where 
Ofc := (1, . . . , 1, a, 1, . . . , 1) G U{N)"^, with a e U{N) in the k-th. place. Then by formula dH), we 
have 

dM = • ^ J] y/Q'{xi)dxi ■ ■ ■ dxNdOi ■ ■ ■ dONdU, 

i 

where 

= I Ai<j,fc (Fij- fc A VFij-fc)|. (5) 

To calculate let us denote by Bi j ^, i ^ j, the derivative of {exp{tEij))i^M (note that since 
Eij lies only in the complexified Lie algebra of U{N)"^, we have (exp(t£'jj))fcM ^ A^mlA, C), 
but this is not important for our considerations). Then equation ([5]) takes the form 

4> = I ^iytj,k Aj,fc|- 

Now (j) can be easily calculated. To do so, we note that for a given i, j, the transformation 
{exp{tEij))k changes only the entries a^j of M. On these entries, it acts by 
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This means that for each i, j, \ Ak Bij±\ = \Jij\, where 



Jij := det 



/ 


^3 





























m 


m 









z^- '■■ 







=n'i 
















s=l 




















V 



















This implies that 

i^j i<j 



Zj Zj I , 



as desired. 



3.3 The probability measure with potential function 
on generalized normal matrices 

Let : C — > M be a potential function. The probability measure on Nm{^,D) corresponding 
to this function is defined similarly to the case of usual normal matrices: 

P^(M)dM = Z-^ exp(-iVtrP^(Mi • • • M„,))dM, 

M G Mmi^, D), where Mi are the blocks of Af . Thus in terms of eigenvalues 



N 



X JJ ^/Q'ixi) JJ \zi - Zjpdxi • • •dxArd6'i • • -d^ArdC/. 
i i<j 

Example 1. Let us calculate the potential function corresponding to the quadratic potential 
Tr (MMt). We have 

Tr(MM^) = XI 141^ = ^{xj+ai) = N^Oi + m^Xj. 

Thus if we choose ai so that Yli «i = (this can be done in a unique way), then 

Tr(MM^) = m^Xj, 
3 

so the corresponding potential function is VF(z) = mQ~^{\z\^) (the function Q is invertible on 
the interval [—a, 00), where a = minaj). 



4 Equilibrium measure 

4.1 Some basic facts about equilibrium measure 

Let Z) be a compact subset of the complex plane C, and l^(z) a potential function (a continuous 
function on D). Denote by ^A{D) the set of the Borel probability measures cj on D without 
point masses, and define the energy of a to be 

la = W{z)da{z)+ / / \og\z - w\~^da{z)da{w). 
Jd Jd Jd 
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An equilibrium measure for 1^ on Z) is a measure a G M{D) such that 
la = inf Jo. 

Theorem 3. The equilibrium measure a exists and is unique. It satisfies the equation 

W{z) - 2 1 \og\z- w\da{w) = C, (6) 
Jd 

where C is a constant, almost everywhere with respect to a. 
The proof of this theorem can be found in [3] . 

Note that equation ([6]) does not have to hold outside the support of a. 

Note also that if a is absolutely continuous with respect to the Lebesgue measure near 
a point zo in the interior of D, and da = g{z)d'^z, where g is continuous near zq and g{zo) > 0, 
then AW = Ang near zq. This clearly cannot happen at points where ISW < 0. In particular, 
if ISW < everywhere, then do" tends to be concentrated on the boundary of D. 

4.2 Asymptotic eigenvalue distribution in the normal matrix model 

In Section [21 we defined a measure 

PN{M)dM = Jn{zi, ZN)d'^zi ■ ■ ■ d'^ZNdU. 
by formula ([2]). We are interested in the behavior of this measure when N oo. Let 5z = 

N 

jj ^ 6z he the measure on D corresponding to the points zj. Then 

- log(Z^ Jjv(zi, . . . , zn)) = (/ WiOddziO + 11^^^ log 1^ - Cr'd6z{C)dSM 

This shows that the leading contribution to the integral with respect to the measure P/v(M)dM 
comes from configurations of eigenvalues zi, . . . ,zn for which the expression in parentheses in 
the last equation is minimized. This means that in the limit ^ oo, we should expect the 
measures 5z for optimal configurations to converge to the equilibrium measure with potential 
function W. 

This indeed turns out to be the shown by the following theorem, proved in [3]. 

Theorem 4. Let the k-point correlation function be 




Then the measure 

RP{{z^)Ufld'z. 

1=1 

on converges weakly to da^^, where da is the equilibrium measure on D, corresponding to 
the potential function W. 

In particular, if = 1, it means that the eigenvalue distribution tends to the equilibrium 
measure in L> as A ^ oo. 
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4.3 Asymptotic eigenvalue distribution 

in the generalized normal matrix model 

As we have seen above, the eigenvalue distribution in the generalized normal matrix model is 



In the limit N ^ oo the second term becomes unimportant compared to the first one, which 
implies that Theorem U] is valid for the generalized normal matrix model. Thus in the limit 
— > oo, the usual and the generalized normal matrix models (with the same potential) are 
equivalent. 

5 Reconstruction of the boundary of the domain 

In previous sections, we showed that in the normal matrix model and the generalized normal 
matrix model, when N ^ oo, the eigenvalue distribution converges to an equilibrium measure 
on D corresponding to some potential function W. In this section, we will try to find this 
measure explicitly in some special cases. 

More specifically, we will consider the case when AW > 0. In this case, if the region D 
is sufficiently large, it turns out that the equilibrium measure is often absolutely continuous 
with respect to Lebesgue measure, and equals da = {'iTT)~^XE^Wd'^ z, where E is a region 
contained in D (the region of eigenvalues), and xe is the characteristic function of E. More 
precisely, it follows from Proposition 3.4 in [3] that if there exists a region E C D such that 
da satisfies equation ([6]) in E, and the left hand side of this equation is > C on D \ E, then 
do" is the equilibrium measure on D for the potential function W. Moreover, note that if E 
works for some D then it works for any smaller D' such that E C D' C D. So, in a sense, E 
is independent of D. (Here we refer the reader to [5], section 4, where there is a much more 
detailed and precise treatment of equilibrium measures, without the assumption A$ > 0). 

Thus let us assume that E exists, and consider the problem of finding it explicitly given the 
potential W. 

5.1 The reconstruction problem 

We will consider the case when D = D{R) is the disk of radius R centered at the origin, and 



where <I> is a function of one variable continuous on [0, oo) and twice continuously differentiable 
on (0, oo), and P a complex polynomial. We assume that (s<&'(s))' is positive, integrable near 
zero, and satisfies the boundary condition lim^-to s*I*'(s) = 0. Computing the Laplacian of W, 
we get (taking into account that A = 495): 



Pjv(M)dM = JN,m{zi, ZN)d^Zi ■ ■ ■ d^ZNdU, 



where 




W{z) = <^{zz) - P{z) - P{z) 



g{s) := {A7t)-^AW = ^"^($'(5) + s$"(s)) 



where s = zz. Define the measure da = gd'^z. 
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Suppose that the region E exists, and contains the origin. In this case, differentiating equa- 
tion with respect to z, we have inside E: 

-z^'iz-z) - P'iz) = [ '-^d'w. (7) 
Je z-w 

On the other hand, inside the disk D, the function 
Wo{z) := 2 / g{ww) log \z — w\d'^w 



JD 

satisfies the equation AWq = 47r(7, and is rotationahy invariant, so 

Wo{z) = ^{zz) + C\ 
where C is a constant. Hence, differentiating, we get, inside D: 

-z^\zz) = [ '-^d^w. (8) 
Jd z-w 

Thus, subtracting ([7]) from ([8]), we obtain inside E: 

P'iz) = [ ^d^«;. (9) 
Jd\e z-w 

Let I{s) = TT J^g{t)dt = s^'{s). Then dl{zz) = -nzg{zz)dz. Thus, using Green's formula, we 
get from ([9|): 

T>i( \ ^ f I{ww) 

P{z) = — / — -dw, 

27rz Jao-dE w{z - w) 

where the boundaries are oriented counterclockwise. The integral over the boundary of D is 
zero by Cauchy's formula, so we are left with the equation 

P'{z) = / 4^d«;. 



27ri Jqe w{w — z) 

This equation appeared first in the theory of Hele-Shaw flows on curved surfaces in and 
it can be solved explicitly by the method of singular points developed in the same paper. Let 
us recall this method. 

5.2 The singular point method 

Define the Cauchy transform oi E with respect to the measure dcr by 

h,{z)=j z^E. 
Jd\e z-w 

This is a holomorphic function of z which (as we have just seen) is independent of the radius R 
of D. As we have seen, it is also given by the contour integral 

hE{z) = — / '-^dw, 10 

2-Ki Jqe w{w - z) 

and in our case we have hE{z) = P'{z). 

Let / : D{1) ^ C \ S be a conformal map, such that /(O) = oo, and (l//)'(0) = a G M+ 
(such a map is unique). 
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Lemma 2. The function 



m=l{f{C)f{Q))-f{Q)hE{f{Q)) 

continues analytically from the unit circle to a holomorphic function outside the unit disk. 

Proof. By the Cauchy formula, we have 

If hE(w) 

hE{z) = / dw, for any z £ E. 

2m Jqe w- z 

So by formula (jlOp. we have 
1 f I{ww) /w — hE{w) 



-dw = 0, for any z £ E. 

2'Ki Jqe w — z 

It follows that the function I{zz)/ z — hE{z), defined along dE, can be analytically continued to 
a holomorphic function outside E, which vanishes at infinity. This implies the lemma. ■ 

Similarly to [3], this lemma implies the following theorem. 



Theorem 5. The function He is rational if and only if the function 9{(^) = I{fiC)f{^/C)) ^s. 
Moreover, the number of poles of 9 is twice of the number of poles of zhE{z). More specifically, 
if (q and 1/Co ctre poles of 9 of order m, then zq = /(Co) ^-5 a pole of order m for hE{z), and 
vice versa. 

Thus, if is a rational function, then 9 can be determined from h up to finitely many 
parameters. 

After this, / can be reconstructed from 9 using the Cauchy formula. For this, note that the 
function / is invertible, since I' = g > 0. Also, 6 takes nonnegative real values on the unit circle. 
Thus, we have 



/(C)/(i/C) = /"'(^(C)). 

Taking the logarithm of both sides, we obtain 



iog(C/(C)) + log (C-V(i/C)) = log (0(C)). 

Thus we have 



2vri 7|„|=i n - C 

1 f logI~^{9{u)) ^ ' 
a = exp I —737 / dn 

|?i|=i , 



The unknown parameters of 9 can now be determined from the cancellation of poles in 
Theorem [5l similarly to the procedure described in [lOj . We note that the knowledge of the 
function He is not sufficient to determine E (for example if S is a disk of any radius centered 
at the origin then He = 0). To determine the parameters completely, we must also use the 
information on the area of E: 



Je 27rz f{u) 
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5.3 The polynomial case 

In particular, in our case, 

hEiz) = P'{z) =ai + a2Z-\ h adz'^'^, 



which imphes that 0(C) = ^jC'^ ^ and bj 

j=-d 

So we get 



-r 



( 



f{C)=aC'exp 



1 

27ri 



log /-I Z bju^ 



u\=l 



u-C 



-du 



a = exp 



1 
47ri 



logr 



-du 



u\=l 



U 



V 



(11) 



Finally, note that if the coefficients of the polynomial P are small enough, then all our 
assumptions are satisfied: the region E exists (in fact, it is close to a disk), and contains the 
origin. Also, in this case the left hand side of equation ([6]) is > C, which implies that the 
equilibrium measure in this case (and hence, the asymptotic eigenvalue distribution) is the 
measure da in the region E. 

Example 2. Consider Example [TJ the generalized normal matrix model with the density 
exp (-/3tr(MMt - P{M) - P{M)^). As we showed, in this case ^{s) = mf3Q (s). So a short 
computation shows that 



TTm~^(3~^g{Q{x)) 



Q{x) 



1 



This implies that g > 0, i.e. our analysis applies in this case. 



5.4 Some explicit solutions 

Consider the case ^{s) = Cs^, C,b > 0. For example, in the generalized normal matrix model 
with at = and potential term as in Example [U one has $(s) = ms^/^, which is a special case 
of the above. 

We have g{s) = 7r~^Cb'^s^~^ , so our analysis applies (note that if 6 < 1 then g is singular 
at zero, but the singularity is integrable and thus nothing really changes in our considerations), 
and I{s) = Cbs^ . Thus the integral in pT]) can be computed explicitly (by factoring 9), and the 
formula for the conformal map / simplifies as follows: 

d 

i=i 

The parameters a > and C,j are determined from the singularity conditions and the area 
condition. 

Consider for simplicity the example d = 1. In this case we have 



hE{z) = K, 
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and we can assume without loss of generality that iiT G R. Then 
/(C) = (aC)-^(l + /3C)'/', /3GM, 

and 

e{C) = Cba~^\l + f3C){l + PC^)- 
The residue of 9 at zero is thus Ch(5a~'^^. Thus the singularity condition says 

Ch(5a^-'^^ = K. 
The area condition is 

l = C76a-2''(l + /?2(i_6-i)). 

Thus we find (3 = KC~^b~^a'^''~^ , and the equation for a has the form 

Cba-^' + C-'b-'K^a''-'{l - 6"^) = 1. 

Remark 2. This example shows that to explicitly solve the generalized (as opposed to the 
usual) normal matrix model in the N ^ oo limit with the quadratic (Gaussian) potential, one 
really needs the technique explained in Section 5 of this paper, and the techniques of [3] are not 
sufficient. 
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